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ABSTRACT 

We present a pilot study of the clustering properties of Lya absorbers with respect to known galax- 
ies based on 112 Lya absorbers and 482 galaxies identified at z < 0.5 along the sightline toward 
PKS0405— 123. The principal goal is to determine the origin of Lya absorbers based on their cross- 
correlation amplitude with known galaxies and investigate a possible N(H I) dependence of the cross- 
correlation function. The main results of our study are as follows. (1) The cross-correlation function £ ga 
measured using only emission-line dominated galaxies and Lya absorbers of log A(H I) > 14 shows a 
comparable strength to the galaxy auto-correlation function £ gg on co-moving, projection distance scales 
< lh Mpc, while there remains a lack of cross-correlation signal when using only absorption- line 
dominated galaxies. This signifies a morphology-dependent £ gQ and indicates that strong absorbers of 
log N(H I) > 14 and emission-line galaxies reside in the same halo population. (2) A maximum-likelihood 
analysis shows that Lya absorbers of log A(H I) < 13.6 are consistent with being more randomly dis- 
tributed with respect to known galaxies. Finally, (3) we find based on this single sightline that the 
amplitude of £ ffa does not depend sensitively on A(H I) for strong absorbers of log N(R I) > 13.6. 

Subject headings: galaxies: evolution — quasars: absorption lines 
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1. INTRODUCTION 

The forest of Lya absorption line systems observed in 
the spectra of background QSOs offers a sensitive probe 
of the tenuous, large-scale baryonic structure that is oth- 
erwise invisible (e.g. Rauch 1998). Over the last several 
years, various numerical simulations have predicted that 
approximately 40% of the total baryons at redshift z = 
reside in diffuse intergalactic gas of temperature T < 10 5 
K that gives rise to the Lya absorbers (e.g. Dave et al. 
2001). Indeed, Penton, Stocke, & Shull (2002, 2004) argue 
that Lya absorbers of neutral hydrogen column density 
N(H I) < 10 14 - 5 cm" 2 may contain between 20 - 30% of 
the total baryons in the nearby universe based on obser- 
vations of the local Lya forest and a simple assumption 
of the cloud geometry. Furthermore, their comparison be- 
tween galaxies and absorbers along common lines of sight 
implies that roughly 20% of low-redshift Lya absorbers 
originate in regions where no luminous galaxies are found 
within Ihr 1 Mpc radius. They conclude that « 5% of all 
baryons are in voids. 

Whether or not Lya absorbers trace the typical galaxy 
population bears directly on the efforts to locate the miss- 
ing baryons in the present-day universe (Persic & Salucci 
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1992; Fukugita, Hogan, & Peebles 1998), and to apply 
known statistical properties of the Lya absorbers for con- 
straining statistical properties of faint galaxies. This issue 
remains, however, unsettled (Lanzetta et al. 1995; Stocke 
et al. 1995; Chen et al. 1998, 2001; Penton et al. 2002). 
While nearly all galaxies within 180 kpc physical ra- 
dius of the QSO lines of sight have associated Lya ab- 
sorbers of JV(H I) > 10 14 cm -2 that are less than a ve- 
locity separation v — 250 km s -1 away along the sightline 
(Chen et al. 1998, 2001), not all Lya absorbers have a 
galaxy found within 1 hr 1 Mpc physical distance to a lu- 
minosity limit of 0.5 L* (Morris et al. 1993; Tripp et al. 
1998). Despite the discordant interpretations in associat- 
ing Lya absorbers with galaxies, the Lya absorbers (par- 
ticularly those of A(H I) > 10 14 cm~ 2 ) are consistently 
found to be correlated with known galaxies (Morris et al. 
1993; Bowen, Pettini, & Blades 2002) but at a weaker clus- 
tering strength than what is measured between field galax- 
ies. The weak cross-correlation amplitude between galax- 
ies and absorbers has led to the conclusion that the Lya 
absorbers trace large-scale filamentary structures, rather 
than individual galactic halos (e.g. Grogin & Geller 1998; 
Tripp et al. 1998). 

However, galaxies are known to be a biased tracer of 
the underlying dark matter halos and the bias depends on 
morphology and luminosity (see e.g. Davis & Geller 1976; 
Loveday et al. 1995; Madgwick et al. 2003). The weaker 
cross-correlation amplitude observed between galaxies and 
Lya absorbers, when compared with galaxy auto-correlation 
studies, can possibly be explained by a difference in galaxy 
bias. In this Letter, we present a pilot study of the galaxy- 
Lya absorber cross-correlation function based on the Lya 
absorbers and galaxies identified along the sightline to- 
ward PKS0405— 123. The primary objectives of the study 
are (1) to compare the galaxy-absorber cross-correlation 
function with galaxy auto-correlation function for galax- 
ies of different properties, and (2) to obtain a quantitative 
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Fig. 1. — Left: The top panel shows the redshift distribution of all 112 Lya absorbers identified along the sightline toward PKS0405— 123. 
The detection limit of Lya absorption lines is uniform at W(H I) = 13.3 across the entire redshift range over z = 0.002 — 0.423 for absorbers 
with Doppler parameter b < 40 km s —1 , and better than 7V(H I) = 13.1 over z = 0.02 — 0.243. Lya absorbers at z > 0.423 are found in 
the FOS spectra retrieved from the HST data archive. The bottom panel shows the cumulative 7V(H I) distribution of the absorbers. Right: 
Redshift distribution of galaxies with R < 20 in the field around PKS0405— 123. The top panel shows all galaxies over the 1600 arcmin 2 
sky area centered at roughly at the background QSO, for which we have obtained robust redshift measurements. The bottom panel shows 
galaxies within 1 h^ 1 Mpc comoving radius in each redshift bin. In both panels, the open histogram represents galaxies of all spectral types, 
the shaded histogram represents emission-line dominated galaxies, and the solid curve shows the selection function of our redshift survey, 
which is 90% complete at R = 20. 



estimate of the N(H I) threshold below which the Lya ab- 
sorbers no longer trace the large-scale galaxy structures. 

2. THE LYa-ABSORBER AND GALAXY SAMPLES 

Exquisite cchcllc spectra of the QSO PKS0405-123 (z qso = 
0.5726) have been obtained by the STIS GTO team (PID = 
7576), using the Space Telescope Imaging Spectrograph 
(STIS) on board the Hubble Space Telescope (HST) with 
a 0.2" x 0.06" slit and the E140M grating (R = 45, 800 
or 6.7 km s _1 ) for a total exposure time of 27,208 s. The 
final stacked spectrum covers a spectral range spanning 
from 1140 to 1730 A and has on average S/N « 7 per 
resolution element over the entire spectral region. Ad- 
ditional spectra obtained using the Faint Object Spectro- 
graph (FOS) have been included for finding strong absorp- 
tion lines log N(H I) > 14 outside of the STIS/E140M 
spectral range at 0.423 < z^ ya < 0.557. Detailed data re- 
duction and absorption line measurements are presented 
in Williger et al. (2005). We have identified 112 Lya ab- 
sorption lines with log TV (HI) = 12.5—16.5 at greater 
than 4-cr significance level over 0.01 < z-^ ya < 0.54. The 
redshift distribution and cumulative N(R I) distribution 
functions of the Lya absorber sample along this sightline 
are presented in the left panels of Figure 1. 

We have carried out a redshift survey of galaxies in the 
field around PKS0405-123, using the WFCCD multi-slit 
spectrograph mounted on the 2.5 m du Pont telescope 
at the Las Campanas Observatory. The primary goal of 
the project is to study the large-scale correlation between 
galaxies and QSO absorption- line systems at low redshift. 
We therefore aim to obtain a complete, magnitude-limited 



survey of galaxies over a field size that is comparable to 
the characteristic correlation length of field galaxies at z ~ 
0.2. Detailed descriptions of the galaxy survey and reduc- 
tion of the spectroscopic data are presented in Prochaska 
et al. (2005, in preparation). In summary, we have ob- 
served 535 galaxies of i?-band magnitude R < 20 over a 
1600-arcmin 2 non-contiguous region centered on the back- 
ground QSO. Robust redshift measurements are available 
for 482 of these galaxies. Our survey is sensitive to galax- 
ies of 0.4 L* at z = 0.2 within 5 h~ x Mpc comoving radius 
from the QSO lines of sight. t 

The redshift distribution of all 482 galaxies identified 
in the field around PKS0405— 123 is presented in the top- 
right panel of Figure 1 . Redshifts are determined using 
multiple spectral line features identified in the data. Typi- 
cal redshift measurement errors are Az = ±0.0002 or ~ 40 
km s _1 . We also perform a x 2 analysis, which compares 
the galaxy spectra with a model established from a linear 
combination of four principal components from the Sloan 
Digital Sky Survey galaxy sample (Schlcgcl et al. 2005, 
in preparation). This procedure allows us to character- 
ize each galaxy based on its emission- and absorption-line 
properties. We found that 284 of all galaxies in our sur- 
vey have emission-line dominated spectral features. The 
redshift distribution of 61 galaxies found within a comov- 
ing radius of p = 1 hT 1 Mpc from the QSO line of sight 
is presented in the bottom-right panel of Figure 1 . The 
shaded histogram represents emission-line galaxies, while 

t Throughout the paper, we adopt a A cosmology, SJm = 0.3 
and S7a = 0.7, with a dimensionless Hubble constant h = 
ff /(100 km s- 1 Mpc" 1 ). 
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the open histogram represents all galaxies identified in the 
same volume. The solid curve in each panel shows the 
selection function of our survey. 

3. ANALYSIS 

3.1. The Galaxy-Absorber Cross-correlation Function 

In the formalism developed by Mo & White (1996), 
the halo-mass cross-correlation function £/im( r ) is related 
to the mass auto-correlation function £ TOm (r) by a mass- 
dependent bias b(M), as £/j m (r) = b(M) £ TOm (r). Compar- 
ing the amplitudes of the galaxy-absorber cross-correlation 
function £ ga and the galaxy auto-correlation function £ gg 
therefore yields an estimate of the absorber halo mass, if 
the mean halo mass of the galaxies is known. This ap- 
proach, which has been applied to study the origin of 
intcrmcdiate-redshift Mgll absorbers (Bouche, Murphy, 
& Peroux 2004), has an important advantage in that the 
complicated selection function of a galaxy survey impacts 
both £ ga and £ gg in the same way and is therefore cancelled 
out in the calculation. 

In principle, the correlation function is multi-variate. It 
depends on the velocity separation between galaxy-absorber 
or galaxy-galaxy pairs v along the line of sight, the comov- 
ing projected distance of a pair p, and the underlying halo 
mass of the source, for which we use the spectral properties 
of the galaxies as a proxy. In agreement with previous find- 
ings of Lanzetta, Webb, & Barcons (1997), who considered 
£ go as a function of v and p, we find that ^ ga (r>) measured 
from the sample of 112 Lya absorbers and 482 galaxies 
identified along the sightline toward PKS0405— 123 weak- 
ens toward larger p and shows no signal beyond p = 1 h~ x 
Mpc. Here we focus on a 1-D correlation function mea- 
sured versus v, and consider only pairs at p < 1 h~ x Mpc. 
We compare £ ga (v) measured for Lya absorbers of differ- 
ent N(H I) with £ gg {v) measured separately for all galax- 
ies, absorption-line dominated galaxies, and emission-line 
dominated ones (Figure 2). Error bars in each velocity 
separation bin are estimated assuming Poisson statistics. 

3.2. Clustering of Lya Absorbers with Galaxies of 
Different Spectral Morphology 

Three distinct features are qualitatively apparent in the 
left panels of Figure 2. First, a strong cross-correlation 
signal in £ ga (v) for Lya absorbers of log N(R I) > 13.5 is 
present over velocity separation v < 250 km s _1 . The sig- 
nal has > 4 — a level of significance when compared with 
random galaxy and absorber pairs found at large veloc- 
ity separation along the sightline, confirming that these 
absorbers are not randomly distributed with respect to 
galaxies. Second, this cross-correlation signal in £ sa (i>) 
becomes weaker for absorbers of lower iV(H I), suggesting 
a correlation between 7V(H I) and the mass scale of the un- 
derlying dark matter halos traced by these absorbers. In 
addition, the absence of a cross-correlation signal for Lya 
absorbers of log N(R I) < 13.5 indicates that these weak 
absorbers do not trace the spatial distribution of galaxies. 
Third, we find that when considering all galaxies together 
the amplitude in £ gg is more than 50% stronger than the 
amplitude in £ ga , even for the strong absorbers in the top 
panel. This is consistent with previous findings, which 
have been interpreted as Lya absorbers tracing large-scale 



filamentary structures rather than individual galactic ha- 
los (e.g. Grogin & Geller 1998; Tripp et al. 1998). 

To examine the nature of Lya absorbers, we repeat 
the analysis including galaxies of different spectral prop- 
erties. We divide the galaxies into two subsamples de- 
pending on whether their spectra are dominated by ab- 
sorption or emission features. Our choice is guided by 
the fact that absorption-line dominated galaxies are more 
strongly clustered than emission-line dominated galaxies 
(e.g. Madgwick et al. 2003). The results of the correlation 
analysis including only absorption-line and only emission- 
line dominated galaxies are presented respectively in the 
middle and right panels of Figure 2. We find that £ ga .{v) 
measured using only absorption-line dominated galaxies is 
consistent with zero on all velocity separation scales for the 
entire Lya absorber sample, contrary to the enhanced £ gg 
measured for absorption-line dominated galaxies only. In 
addition, we find that £, ga (v) = 6.2 ± 1.9 for Lya absorbers 
of log N(K I) > 14 and emission-line galaxies are consis- 
tent with t; gg (v) = 6.6 ± 0.8 for emission-line galaxies to 
within 1-sigma uncertainties, indicating that these strong 
Lya absorbers and emission-line dominated galaxies found 
in our redshift survey trace the same halo population. The 
weaker £ ga (v) at lower N(R I) observed in the right pan- 
els also suggests that weaker absorbers may originate in 
still lower mass halos that are loosely clustered with these 
known galaxies. 

3.3. N(H l)-Dependent Clustering Strength of Lya 
Absorbers with Galaxies 

Much of the past dispute regarding the nature of Lya 
absorbers originates in the different iV(H I) regimes con- 
sidered by different authors. It appears that while most 
strong absorbers of N(H I) > 14 can be explained by the 
known galaxy population (e.g. Chen et al. 1998, 2001), 
the majority of weaker absorbers remain as a poor tracer 
of galaxies (e.g. Penton et al. 2002, 2004). This is also 
observed in the declining cross-correlation amplitude with 
iV(H I) as presented in Figure 2. We perform a maximum 
likelihood analysis to obtain a quantitative estimate of the 
N(R I) threshold below which the Lya absorbers no longer 
trace known galaxies. 

We first formulate the cross-correlation function as 
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if N(Ll I) > N and £, ga (v) = otherwise. We model the 
velocity dispersion function of the galaxy and absorber 
pairs using a gaussian of FWHM vq. According to the 
model described in Equation (1), we expect that the cross- 
correlation amplitude is scaled with 7V(H I) for absorbers 
of N(R I) > N following a power-law function, and van- 
ishes for weaker absorbers. For each Lya absorber i of 
N l (H I) at z\ jyal the probability of finding a galaxy j 

within the volume defined by z{ jya ± Az and pj < 1 h^ 1 
Mpc is then written as 
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where Vij/c = [(1 + z gal ) 2 - (1 + z LyQ ) 2 ]/[(l + z gal ) 2 + 
(1 + ZLy Q ) 2 ] and S(z) is the redshift selection function 
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Fig. 2. — Galaxy-Ly« absorber cross-correlation function ^ ga (v) versus velocity separation v of each pair along the line of sight in the top 
three panels for Lya absorbers of different N(R I), and galaxy-galaxy auto-correlation function !;gg(v) versus v in the bottom panel. We note 
that the limits on the y-axis for weak absorbers has been rescaled for a better visibility of the cross-correlation signal at v < 500 km s . 
For £ ga , we consider all 61 galaxies found within p = 1 h Mpc co-moving radius to the sightline in the left panels, and 15 absorption-line 
dominated and 46 emission-line galaxies at p < 1 h^ 1 Mpc in the middle and right panels, respectively. For £ g3 , we consider all pairs of 
different spectral type with a co-moving impact separation < 1 h^ 1 Mpc. 



of our galaxy survey within the 1 h~ x Mpc comoving ra- 
dius as presented in the bottom-right panel of Figure 1. 
Considering all the galaxy and Lya absorber pairs as a 
whole, we can formulate the likelihood function of observ- 
ing to galaxies at p < 1 h^ 1 Mpc from n Lya absorbers of 

12.5 < A(H I) < 17 as the following, 

n m 

C(A,N ,(3,v ) = l[l[P i (4j. (3) 

1=1.7=1 

Finally, maximizing C with respect to A, N , (i, and v 
leads to best-fit values of A = 11 ± 3, log N = 13.6 ±0.1, 
(3 = 0.0 ±0.1, and vq = 102 ± 9 km s , where the error 
of each parameter represents the 68% 1-cr, 1-parameter 
uncertainty. 

Adopting a power-law model (A(H Y)/Nq)P with log N = 

13.6 to replace the step function in Equation (1) yields a 
best-fit (3 = 0.20 ± 0.02. The non-zero value reflects a 
significant difference in £ sa between strong and weak ab- 
sorbers. We note, however, that the best-fit power-law 
model substantially underestimates £ 9a at log A(H I) > 
14, when extrapolating from log N(H I) = 12.5 — 13, and 
is therefore not likely to represent the intrinsic shape of 
£, ga - Adopting a gamma-function form following 
exp(-jV /iV(H I)) x (N /N(R I)) 13 yields logA = 13.5 ± 
0.1 and P = 0.0 ±0.1. 

4. DISCUSSION AND CONCLUSIONS 

Our pilot study of the galaxy-Lya absorber cross-correlation 
function at z < 0.5 using galaxies and absorbers identi- 
fied along the sightline toward PKS0405— 123 has offered 
important insights for understanding the nature of Lya 
absorbers over a wide range of A(H I): 

First, the cross-correlation analysis not only confirms 
that Lya absorbers (particularly those of log JV(H I) > 
14) arc not randomly distributed with respect to galax- 
ies, but also shows that emission-line dominated galax- 
ies and these strong absorbers have comparable cluster- 



ing strength. The difference in the auto-correlation am- 
plitudes observed for absorption- and emission-line domi- 
nated galaxies confirms that absorption-line galaxies trace 
higher-mass halos with stronger clustering strength than 
emission-line galaxies. We therefore conclude that Lya 
absorbers of log Af(H I) > 14 and emission-line dominated 
galaxies found in our redshift survey do indeed trace the 
same halo population. A quantitative estimate of the halo 
mass scale traced by the Lya absorbers will require a de- 
tailed clustering analysis of these emission-line galaxies 
over a large volume. 

Second, we observe a likely transition in the cluster- 
ing strength of Lya absorbers with respect to galaxies at 
log A^H I) = 13.6, above which the amplitude of £, ga {v) 
is comparable to the emission-line galaxy auto-correlation 
function and below which ^ ga (v) quickly declines to within 
2.5 a of zero. The diminishing cross-correlation signal be- 
tween these weak absorbers and known galaxies suggests 
that they are a poor tracer of galactic structures. A de- 
tailed comparison with numerical simulations will help us 
to constrain the amount of baryons contained in these low- 
TV (H I) Lya absorbers. In addition, we do not detect a 
strong dependence of £ 9a on A(H I) at log A(H I) > 13.6. 
It would be interesting to see whether £ ga remains insen- 
sitive to A(H I) with a larger galaxy and absorber pair 
sample over multiple lines of sight. 

With a larger sample of galaxies and absorbers over mul- 
tiple lines of sight, we hope to obtain a better constraint 
for the galaxy-Lya absorber cross-correlation function in 
the near future. In addition, we will extend this analysis 
to OVI absorbers identified along these sightlines, in or- 
der to locate the warm-hot intergalactic medium probed 
by these absorbers (e.g. Mulchaey et al. 1996; Tripp, Sav- 
age, & Jenkins 2001). 
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